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Abstract 

Extensive Monte Carlo simulations were carried out to investigate the nature of the ordering 
transition of a model of adsorbed self-assembled rigid rods on the bonds of a square lattice [Tavares 
et. al., Phys. Rev E 79, 021505 (2009)]. The polydisperse rods undergo a continuous ordering 
transition that is found to be in the two-dimensional Ising universality class, as in models where 
the rods are monodisperse. This finding is in sharp contrast with the recent claim that equilibrium 
polydispersity changes the nature of the phase transition in this class of models [Lopez et. al., 
Phys. Rev E 80, 040105(R)(2009)]. 



PACS numbers: 64.60Cn, 61.20.Gy 



I. INTRODUCTION 



It has been shown, recently, that pure hard-rod models in two-dimensions (2D) exhibit 
discrete orientational order without translational order 111, driven by a mechanism resem- 
bling that proposed by Onsager for the nematic transition of rods in three dimensions |2j. 
Specifically, it was proved that a system of rods on the square lattice, with hard-core exclu- 
sion and length distribution between 2 and n, exhibits discrete orientational long-range order 
for suitable fugacities and large n. This may seem surprising as the nature of the transition 
of monodisperse freely rotating rods in 2D remains subtle, as it appears to depend on the 
details of the particle interactions [35]. 

Simple 2D restricted orientation models are relevant to describe the sub-monolayer regime 



of linear molecules adsorbed on crystalline substrates 6j and, even without polydispersity, 
rigid rod (RR) models were shown to exhibit a number of interesting features^- 9] . It was 
found that the ordered phase is stable for sufficiently large aspect ratios [7, 9] and that the 
transition on the square lattice is 2D Ising [81]. 

Polydisperse restricted orientation RR models arc generalizations of the Zwanzig 



model 



101 ]. and provide a useful starting point for understanding the effects of polydispersity 



on the phase behavior of RRs 11|. The description of self-assembled rods has to consider not 
only the effects of polydispersity but also the polymerization process. In this context, we 
proposed a model of self-assembled RR (SARR), composed of monomers with two bonding 
sites that polymerize reversibly into polydisperse chains In the (lattice) model a site 
can be either occupied or unoccupied and each occupied site has a spin variable. On the 
square lattice, the spins take two values representing the discretised set of orientations of 
the bonding sites that coincide with the lattice bonds. The interaction between two spins 
depends not only on their relative orientations but also on their orientations relative to the 
lattice bond connecting the monomers. We used a simple theory to investigate the interplay 
between self-assembly and ordering over the full range of temperature and density. The 
results revealed that the continuous ordering transition is predicted semiquantitatively by 
the theory 12j . The universality class of this transition was not investigated; ordering of 
SARRs was assumed to be that of monodisperse rigid rods, which was found to be 2D Ising 
on this lattice js]. 

In fact, the transition of polydisperse RRs, on the square lattice, was investigated for a 



2 



vertex model that allows configurations promoting the polymerization of rods, in such a way 
that it is equivalent to the hard square model on the diagonal lattice. In polymer language, 
the ordered phase is stable when the average polymer length is long or its density is high. 
Calculations of the order-parameter using a variant of the density matrix renormalization 
group exhibit clear 2D Ising exponents (/3 = 0.125) at all densities[3]. However, in 2009, 
Lopez et. al [l^] carried out Monte Carlo (MC) simulations to investigate the critical 
behavior of the SARR model and concluded that self-assembly affects the nature of the 
transition, claiming it to be in the q=l Potts class (random percolation), rather than in the 
2 D Ising ^ Potts). This is at odds with exact results)] that map the polydisperse RR 
model, with n = oo, to the 2D Ising model, as well as with the results of the vertex model 



131 ] referred to above. 



Apart from its fundamental interest, self-assembly is a very active field of research, driven 
by the goal of designing new functional materials, inspired by biological processes where it 
is used routinely to construct robust supramolecular structures. In this context, the effect 
of polydispersity on the nature of the ordering transition of a given model is an important 
open question. 

In the following, we report the results of a systematic investigation of the criticality of 
the SARR model over a wide range of temperatures, corresponding to critical densities that 
decrease from 1.0 (full lattice) to 0.1. We note that the full lattice SARR model may be 
mapped to a 2D Ising model, while the zero density SARR model exhibits an equilibrium 
polymerization transition at zero temperature. The results of our simulations provide strong 
evidence that the transition remains in the 2D Ising class at all (finite) densities. However, 
the numerical results also suggest that the scaling region is strongly affected by the density, 
decreasing as the density decreases, in a way that depends both on the scaling variable and 
on the thermodynamic function under investigation. 

This paper is organized as follows: the model and the simulation methods are described 
in Sec. II. The results for 2D Ising criticality of the SARR model are reported in section 
III. In section IV we give additional arguments that support our conclusion: (i) we map the 
full lattice limit (FLL) onto the 2D Ising model, (ii) we consider the zero density limit and 
estimate the crossover line from the zero density 'equilibrium polymerization transition' and 
(iii) we discuss the non-monotonic behavior of the internal energy per particle on the critical 
line. Finally, in section V we summarize our results, and offer an explanation for the q=l 
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Potts behavior observed by Lopez et al 



II. THE MODEL 



•0. 



The model is the two bonding site model, on the square lattice, proposed in Ref. [12j in 



the context of a general framework to understand self-assembly (see [15|, |l6j and references 
therein). A lattice site is either empty or occupied by one monomer with two bonding sites. 
Each monomer, i, adopts one of two orientations, Sj = x or Sj = y, corresponding to the 
alignment of the bonding sites with the lattice directions, x and y. Monomers attract each 
other if their bonding sites overlap, promoting the self-assembly of polydisperse rigid rods. 
The energy of the system may be written as: 

M 

U = -e^2 ^2 |s(ri) •s(r i + a)||s(r i ) -a|, (1) 

i=l a=x,y 

where i labels a lattice site, s(r) denotes the monomer orientation (s = for an empty site); 
x and y are lattice unit vectors, and M is the total number of sites. 

The criticality of this model was investigated in Ref. [l4] where it was found that poly- 
dispersity changes the nature of the ordering transition. In order to check this claim we have 
studied the model over a wide range of thermodynamic parameters, using a multicanonical 
MC method based on a Wang-Landau sampling scheme. We considered systems with sizes 
L x = L y = L, M = L 2 sites, and periodic boundary conditions (PBC). The simulation 



methods were used in previous studies and details may be found there|17H19|. Briefly, in a 
simulation run we fix the the system size, L, and the temperature T; we sample over the 
number of particles < N < M and attempt exclusively MC moves of insertion and deletion 
(with equal probability). In the insertion attempts the orientation of the particle is chosen 
at random. The probability of a configuration, R^, with N particles is: 

P(R N \M, T) oc w(N) exp [-U{R N )/k B T] (2) 

where k B is the Boltzmann constant and the function w(N) is chosen to ensure uniform 
sampling of the density. The probability of a configuration with N particles is: 

P(N\T) oc w(N) J dR^exp [-U(R N )/k B T] = w(N)e- A ^ N ' M ^ kBT , (3) 
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where A(N, M,T) is the Helmholtz free energy. The weight function required for uniform 
sampling of N, in the range [0, M], satisfies: 

W(N) ~ e MNM-T)/k B T + x y (4) 

Clearly, the Helmholtz free energy A(N, M, T) is not known a priori, but ap propriate es- 
timates of w(N) may be obtained |l8j using a Wang-Landau-like method [201] . The multi- 
canonical simulation and the computation of the required observables (energy, order param- 
eters, etc.) are then carried out for < iV < M. In line with previous work, we define the 



order parameter as 



12 



\N — N I 

5 = ^iM, (5) 

where N x and N y are the number of monomers oriented in the directions x and y, respectively. 
The ordering transition, at a given temperature, is located by search ing for pseudo-critical 



values of the chemical potential, p c (L,T). We note that Lopez et al.[14| used the density, 
p = N/M, as the control parameter. In the SARR model at fixed T, the chemical potential 
p is the only external field and plays the role of the temperature T in standard (full lattice) 



Potts simulations 2 lj. We proceed by defining analogues of the Ising response functions, 
related with the second derivatives of ^/k B T (with $ = A — Np, the Grand Potential) with 
respect to the coupling constant K = 1/fcgT, and p: 

1 d 2 [K$(p, M, K)} _fd[u- pp] \ 



k B T 2 V dK 2 V dT 



pl _ k B Td 2 [K<f>(p,M,K)\ _ /dp' 



V dp 2 \ dp 



T,M 



where u = U/M. The quantities c, and p'^ are expected to scale at the critical point as [2 11]: 

c(L,p c (T))~L a /\ (8) 

pl(L,ii c (T)) ~ (9) 

where a, and v are the specific heat and correlation length critical exponents. 

We carried out MC simulations at several temperatures. At each temperature, a range 
of system sizes was considered; up to L = 144 at reduced temperatures, T* = k B T/e < 0.25 
and up to L = 112 at higher temperatures. The results of each simulation are used to 
calculate histograms of the different observables that were then computed in terms of the 
chemical potential. 




FIG. 1: (Color on line) Fourth order Binder cumulant as a function of \x for different system sizes, 
at T* = 0.25 and T* = 0.30. 

III. RESULTS 



Binder cumulant 



We start by computing the fourth-order Binder cumulant 

< 5 4 > 

9a 



(10) 



< 5 2 > 2 ' 

as a function of the chemical potential. In Fig. [I] we plot g±{p) for different system sizes, at 
T* = 0.25 and T* = 0.30. It is clear that the cumulants for different system sizes, L, cross 
at a value of that is very close to the universal critical value for 2D Ising systems, with 



PBC and L x = L y (i.e. g\ ~ 1.168) 22]. This immediately suggests that the criticality of 
the SARR model is in the 2D Ising (q=2 Potts) class (Q2UC), in contrast with the findings 
of Lopez et al. [ijj . Similar results were obtained for all the other temperatures investigated. 



B. Computation of (dp/d^) T 

In Fig. H] we plot the derivative of the density with respect to the chemical potential, 
p' as a function of /i, for different system sizes, at the same temperatures T* = 0.25 and 
T* = 0.30. At values of the chemical potential, /x, close to its critical value, the derivative of 
the density, p' exhibits clear signs of singular behavior, with a peak that increases as the 
system size increases. The size dependence of the peaks is analyzed in Fig. [3j The scaling 
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FIG. 2: (Color on line) Derivative of the density, p, as a function of the chemical potential, /i, 
for different system sizes, at T* = 0.25 and T* = 0.30. A singularity is clearly signalled at both 
temperatures. 

T*= 0.25 T*=0.30 




FIG. 3: System size dependence of the maximum of kBT(dp/dp)T, at T* = 0.25 and T* = 0.30. 
Simulation results are represented by points. Lines are fits to the scaling laws discussed in the text. 

of the peaks of (p') rnax (L, T), with the system size L, is characteristic of the universality 



class of the transition 21]. For SARR on the square lattice we anticipate either q=l Potts 

n n 

(Q1UC) behavior {a/v = —1/2) [23( as reported in Ref. [14j or Q2UC behavior, as found 
for monodisperse rods on the same lattice[8|. In the latter case ajv = and the peak is 



expected to diverge logarithmically 



231 ] . The two scaling laws are: 



U ) max {L) =a + a.L- 1 ' 2 , for Q1UC, 



(p') maX ( L ) = a o + «i In L, for Q2UC. 



(11) 
(12) 
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FIG. 4: (Color on line) System size dependence of the pseudo-critical chemical potentials. See the 
text and the legend for details. Lines are least-squares fits to the MC results. 

In Fig. |3l we plot fits of the two scaling laws to the simulation data. In both cases we 
discarded the data for the smallest system sizes. We found that the simulation results are 
better described by the 2D Ising scaling law, at all temperatures. In fact, equation f lT2|) fits 
the data over a broader range of system sizes (L > 24 at T* = 0.25 and L > 32 at T* = 0.30) 
than does Eq. CD] (L > 56 at both temperatures). 



C. Critical Line 



We start by defining the pseudo-critical chemical potentials at fixed temperature. We 
consider the Binder cumulants, as functions of the chemical potential, and define the pseudo- 
critical chemical potentials, = /J, C (L) (at given T), such that: 

g,(L,^ c L \T)=gl (13) 

We have also used different definitions based on the position of the maxima of the density 
fluctuations and of (<91n < S > /dfj,)?, to check the consistency of the results. In Fig. @]we 
plot the pseudo-critical chemical potentials as functions of 1/L, at T* = 0.25 and T* = 0.30. 
At both temperatures, and L > 32, the chemical potential, fi c (L), computed using Eq. (|T3l) 
is almost independent of system size (horizontal lines in the left and right panels of Fig. HI). 
Estimates of the critical chemical potential, fi c , obtained by extrapolating the MC results, 
are collected in Table fl] Considering the behavior of fi c (L) obtained using Eq. (fTBl . we used 
the results of simulations at this chemical potential, to compute the critical density p c and 
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the critical exponents. Assuming 2D Ising behavior, the critical density is given by 211 ] : 

Pc (L jf ii L \T) = p c (T) + aL-\ (14) 



The results for the critical line T c (p) are plotted in Fig. [5j As expected, the temperature at 



the ordering transition decreases as the density decreases. The critical points calculated in 
earlier work (diamonds 12] and square [l^]) fall on the critical line, within the statistical error 



(open circles). The line of critical points of the SARR model continues beyond the lowest 
density reported in Fig. However, the rapid increase of the average length of the rods 
at these (low) densities and temperatures prevents an efficient simulation of these systems 
with the currently available techniques. 
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FIG. 5: (Co 



3, 



and 



14] 



or on line) Critical line of the SARR model. The diamonds and square are from Refs. 
respectively. 



D. Critical Exponents 

The critical exponents /3/v, and 7/1/ Q, Q were estimated by fitting the MC results to 
the scaling laws, 

\og5{L,pf\T) = a p -^\ogL 1 (15) 

and 

log X (L, pf\ T) = a € + ^ log L, (16) 



where the susceptibility x is defined as: 

X = L 2 [< 5 2 > — < 5 > 2 ] fk B T. (17) 
The critical parameters and effective exponents are collected in Table [J at 10 different 
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0.9338(4) 
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0.123(2) 


1.751(4) 
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oo 


1.000 


-0.85355 


0.125 


1.750 



TABLE I: Results for the critical parameters and effective critical exponents. 

temperatures. Note that the exponents computed for (3/v lie between those corresponding 
to the Q1UC (/3/u = 5/48 ~ 0.104)[l^, Q at low temperatures and those corresponding to 
the Q2UC (J3/v = 1/8) [23] at high temperatures. The results for y/v are closer to those of 
Q2UC (y/u = 7/4) over a wider range of temperature but at the lowest temperatures they 
also approach those of Q1UC (y/u = 43/24 ~ 1.792). 

We note that Eq. (j!5p does not provide a good fit of the simulation results for a wide 
range of system sizes, and thus the results for /3/u should be regarded as effective exponents. 
Consideration of higher-order finite size corrections, of the form, 5 C (L) = L~^' v [ao + bL^], 
is not a feasible, as fits of the simulation results with four adjustable parameters cannot 
discriminate between these two, similar, scaling laws. 

To proceed, we have investigated the finite-size scaling of the Binder cumulant, at T* = 
0.15 and T* = 0.40. In Fig. [6] we plot g^{n) versus L x l v {ii — /i c ), using for the critical 
exponent, u, the values corresponding to the q=l [u = 4/3) and q=2 [u = 1) universality 
classes. At T* = 0.40, the data collapse with v = 1 is excellent, confirming that scaling for 
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the Q2UC is satisfied for all the system sizes (L > 48). However, at the lowest temperature, 
T* = 0.15, the data collapse fails for both universality classes and systems with L < 80. 
Nevertheless, the collapse observed with the Q2UC exponents is marginally better than that 
observed with the Q1UC exponents, suggesting that the SARR criticality is still in the 2D 
Ising class. 
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FIG. 6: (Color on line) Fourth order Binder cumulant as a function of \x for different system sizes, 
at T* = 0.15 and T* = 0.40. 



Finally, we analyzed the scaling behavior of the derivative of the logarithm of the order 
parameter with respect to the chemical potential, at constant temperature. The maxima of 
this quantity are expected to scale with the system size as [l^,[2l|, 24] Q'^ = (d\n5/d/i)^: a ^ ~ 
L x l v . On the basis of this scaling law, we have computed effective values of 1/V by Chi- 
square fitting 25] the simulation results to Q' [V) = aL l l u . The results for \/u are collected 
in table [Til f° r several temperatures, and confirm that the SARR model is in the 2D Ising 
class. It is also clear that as the temperature decreases the effect of the finite system size 
becomes more important, i.e. one requires larger systems to stay on the asymptotic scaling 



region. 
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n 




^max 


X 2 /d.o.f 


1/v 


0.15 
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80 


144 


0.84 


1.22(15) 


0.20 


5 


80 


144 


1.71 


1.13(10) 


0.25 


5 


80 


144 


1.13 


1.04(9) 


0.30 


8 


18 


112 


0.08 


1.05(4) 


0.35 


8 


10 


112 


0.65 


1.02(2) 


0.40 


10 


32 


112 


1.33 


1.00(2) 



TABLE II: Estimates of effective values of 1/v for different temperatures; n is the number of points 
(system sizes) used in the fitting; L m i n , and L max are the minimum and maximum system sizes 
considered, with L m i n chosen to provide statistically acceptable values for the \ 2 merit function 
251 ] . d.o.f. is the number of degrees of freedom in each fitting. 



IV. ANALYSIS OF THE CRITIC ALITY OF THE SARR MODEL 
A. The full lattice limit: 2D Ising 

The results of the previous section suggest clearly that the criticality of the SARR model 
is in the 2D Ising class. In this section we investigate the full lattice limit of the SARR 
model, or full lattice limit (FLL) for short, where we can prove that this is indeed the case. 



One also expects t 
on the critical line 



le criticality to remain unchanged, as long as no other transitions occur 



26] 



We performed a number of simulations using a multi-temperature algorithm proposed by 
Zhang and Ma 27j] and found that the critical temperature is the same as that of the lattice 
gas. This is more than a coincidence as shown below. 

We have mapped the monomer orientations x, y to the Ising spins ±1 and computed the 
total energy of the models by adding the contributions, u p , of square elementary plaquettes, 



u 



where each plaquette consists of a square with four sites enclosing an elementary cell of the 
lattice, where we have taken into account that each pair interaction is counted in two different 
plaquettes. In Table HTT1 we collect the energies for representative plaquette configurations of 
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Plaquette 


+ + 
+ + 


+ + 
+ - 


+ + 


+ - 
- + 


Ising 


Up/ei 


-4 








4 


SARR 


u p /e 


- 2 


-1 


- 1 






TABLE III: Plaquette interactions in the Ising and the FLL of the SARR models on the square 
lattice. 

both models. The mapping between the two models, is then for any plaquette configuration: 

u J" n 7ej = 4< Li /e + 4- (19) 



implying that the FLL limit of the SARR model, on the square lattice, is in the 2D Ising 
universality class. 



B. The zero density limit: Self-assembly 

The SARR model has two independent thermodynamic parameters, the temperature 
and the density of monomers. At high temperatures k^T » e there is little bonding 
and the behavior of the model is similar to that of the lattice gas. At low temperatures, 
however, bonding dominates and the model behaves in a strikingly different way. Rods 
self-assemble and, at a fixed density, the average rod length increases exponentially as the 
temperature decreases. The polydisperse rods undergo an ordering transition at a density 
that is temperature dependent. The transition was calculated in Ref. 12| using a generalized 
mean-field theory of self-assembly, where the polydisperse rods interact through Onsager-like 
excluded volume terms only. 



The critical line is given by 



12|: 



1 



In 



(2-p c )(2 + p c ) 
2p3 



(20) 



and is plotted in Fig. [7J This line is singular in the zero-density limit, where the average 
rod length diverges, signalling the self-assembly or equilibrium polymerization transition at 
zero temperature. An estimate of the crossover line, from the polymerzation transition, is 
obtained from the asymptotic relation between T* and p c , at the singular point: 

1 



T ~ 

c 



In p c 



(21) 
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Using the results for the thermodynamic potentials derived in Ref. 12] one finds that 
the chemical potential at the critical point, p c , is given in terms of the critical temperature 
and density, T c and p c : 

fi c /e = T c * [\n(2pl) + p c /2 - 2 ln(2 - p c )] . (22) 
Using the asymptotic form, Eq. (12 ip. we find for the crossover line: 

fi e /e~-l- (23) 

This line delimits the region where the self-assembly or equilibrium polymerization fluctu- 
ations are large. In Fig. [7] we plot the critical and crossover lines of the SARR model as 
functions of p and p. Note that the crossover line approaches the critical line tangentially 
in the equilibrium polymerization limit, suggesting that the asymptotic scaling region of 
the finite density critical point decreases rapidly as the critical temperature and density 
decrease. 

3.0 
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FIG. 7: (Color on line) Critical lines from MC simulation (points), theory of self-assembly (dashed) 
and mean-field theory (dashed-dotted). Crossover line is the full line. Left panel: p,,T diagram. 
Right panel: p, T diagram. See text for details. 



Finally, the internal energy on the critical line is easily calculated and we find 

U_ 2p c 
Ne 2 + p c ' 



12|: 



(24) 



Note that the energy per particle increases monotonically with the (critical) density. This 
is plotted in Fig. and will be discussed in the next section. 
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C. Intermediate densities 



A look at Table H] reveals that the critical energy per particle, u c =< U/N > c , varies 
non-monotonically with the temperature. In the self-assembly limit, the internal energy 
u* =< U/Ne > is a measure of the average rod length and is predicted to vary monotonically 



on the critical line, as stated above (see Ref. 12J for details). Clearly, a departure from this 



behavior indicates the importance of the attractions between rods, which are short in the 
high-temperature regime. 

Let us assume that there is no bonding, i.e., the temperature is so high that the average 
rod length is of order 1. Then, on average, each monomer interacts with its aligned neighbors, 
the pairs being aligned with the corresponding lattice bonds. The mean-field free energy 
becomes: 

0f= J2 p a (lnp a -l) + (l-p)(ln(l-p)-l)-/3e(p^ + P y, (25) 

a=x,y 

where p x and p y are the densities of particles aligned with the x and y directions, respectively, 
and we have accounted for the entropy of the empty lattice sites. Given that p = p x + p y 
and defining A = p x — p y , the free energy may be written in terms of these variables. The 
critical points are obtained by (i) calculating the field associated with A, (5h = ^£ and (ii) 
setting h = to obtain, implicitly, A(p, T). The critical line is given by: 

T* = Pc . (26) 

and the internal energy on the critical line becomes: 

— = ~ P -- (27) 
iVe 2 y ' 

The internal energy on the critical line is plotted in Fig. [HJ The points are the computer 
simulation results while the two lines are obtained from the self-assembly and the high- 
temperature mean-field theories. The results from self-assembly, Eq. ( 12^1) . are monotonically 
increasing, while those from the high temperature theory, Eq. ( 1271) . are monotonically 
decreasing. 



The low density/temperature behavior is captured by the self-assembly theory 12] while 
the high density/temperature limit is described by the mean-field theory. Similar remarks 
apply to the critical line itself, as shown in Fig. [7J 
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FIG. 8: (Color on line) Internal energy per particle on the critical line. Points: MC simulation, 
dashed: theory of self-assembly, dashed-dotted: mean-field theory. See text for details. 

This analysis suggests that although self-assembly fluctuations become increasingly im- 
portant as the density decreases the nature of the singularity changes at p = only. Nev- 
ertheless, the scaling region decreases rapidly in the low-density/temperature region and 
the true asymptotic behavior may be difficult to observe in simulations of reasonably sized 
systems. 

V. DISCUSSION 

The results reported in the previous sections clearly suggest that the critical behavior of 
the SARR model is 2D Ising. This conclusion is supported by (i) the scaling behavior of 
the Binder cumulant for different system sizes, (ii) the system size dependence of the peaks 
of /^(/i) and (iii) the values of the critical exponent v. This conclusion contrasts with that 
of Lopez et al., and it is important to understand the reasons for this discrepancy. First 
we acknowledge that the values of /3/V are relatively similar for the two universality classes; 
the same may be said of 7/V. Thus the distinction between the two universality classes will 
have to be based on the value of the g±{L) crossing and on the value of v. In the analysis of 
Lopez et al, the use of the density as the control parameter leads to a value of the g^ crossing 
that differs substantially from that of the 2D Ising universality class. We have shown that 
using n as the control parameter leads to a more robust scaling of g^ and to a much better 
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FIG. 9: (Color on line) Fourth order Binder cumulants for different system sizes as a function of 
the scaling densities in the q = 1 and q = 2 universality classes, at T* = 0.30. 

overall Ising scaling. 

Concluding that the SARR model is indeed in the Ising universality class, the question is 
then, how was the value of v « 4/3 observed when using p as the control parameter? Con- 
sider a property Q whose derivative with respect to the control parameter has a maximum 
in the critical region. Such derivative scales in the finite size region as 



iL 1/v (l + bL- 



(28) 



where s represents either the density or the chemical potential. In the scaling region Q' 
and Q' p are related by: 

/f),,\ J, 

(29) 



Q'p ~ 



dp 



InL' 



where we used v = 1 (Q2UC) and the scaling relation given in Eq. ( fT2l) . We suspect that 
the presence of InL in the scaling of Q' p , raises the value of the effective critical exponent v. 
In particular, for the range 60 < L < 120 used by Lopez et al. [ijj] the ratio L/\nL is well 
described by: L/lnL ~ aL 1 ^' , with v' ~ 1.291, close to the value v = 4/3 of the Q1UC. 
This is illustrated in Fig. [9] where it is clear that the finite-size scaling of g<i(p, L) is well 
described using both x = L~ 3 / 4 (p — p c ) and x = (p — p c )L/ InL as the scaling density. To 
conclude, we have shown that the criticality of the SARR model is 2D Ising. Nevertheless, 
as the temperature decreases deviations from the Ising scaling laws increase, and larger 
system sizes are needed to obtain accurate estimates of the critical exponents. This can be 
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understood in terms of the self-assembly fluctuations that occur closer to the critical line 
as the density and temperature decrease. In addition, the use of PBC may enhance this 
finite-size effect, through the percolation of 'periodic' rods. 

analysed the critical behavior of a flexible 



Finally, we note that Milchev and Landau [28] 
self- assembling rod model in 2D. They report a continuous transition in the T, \i space ending 
at a tricritical point, at finite density, and critical exponents on the continuous portion of 
the two-phase boundary in the 2D Ising class. Their model is richer than ours but the nature 
of the continuous portion of the phase boundary is likely to be the same. The connection 
between these two models as well as extensions to 3D will be left for future work. 
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